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Context. Recent observations revealed that the solar atmosphere is highly structured in density, temperature and magnetic field. 
The presence of these gradients may lead to the appearance of currents in the plasma, which in the weakly collisional corona 
can constitute sources of free energy for driving micro-instabilities. Such instabilities are very important since they represent a 
possible source of ion-cyclotron waves which have been conjectured to play a prominent role in coronal heating, but whose solar 
origin remains unclear. 

Aims. Considering a density stratification transverse to the magnetic field, this paper aims at studying the possible occurrence 
of gradient-induced plasma micro-instabilities under typical conditions of coronal holes. 

Methods. Taking into account the WKB (Wentzel-Kramers-Brillouin) approximation, we perform a Fourier plane waves analysis 
using the coUisionless multi-fluid model. By neglecting the electron inertia, this model allows us to take into account ion- 
' cyclotron wave effects that are absent from the one-fluid model of magnetohydrodynamics (MHD). Realistic models of density 
and temperature, as well as a 2-D analytical magnetic-field model, are used to define the background plasma in the open-field 
funnel in a polar coronal hole. The ray-tracing theory is used to compute the ray path of the unstable waves, as well as the 
evolution of their growth rates during the propagation. 

Results. We demonstrate that in typical coronal hole conditions, and when assuming typical transverse density length scales 
taken from radio observations, the current generated by a relative electron-ion drift provides enough free energy for driving the 
mode unstable. This instability results from a coupling between oppositely propagating slow-mode waves. However, the ray- 
tracing computation shows that the unstable waves propagate upward to only a short distance but then are refiected backward. 
k>( ^ The corresponding growth rate increases and decreases intermittently in the upward propagating phase, and the instability 
■ ceases while the wave is propagating downward. 
Conclusions. Drift currents caused by fine density structures in the magnetically open coronal funnels can provide sufficient 
energy for driving plasma micro-instabilities, which constitute a possible source of the ion-cyclotron waves that have been 
invoked for coronal heating. 
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1. Introduction 

New observations of coronal structures made by the high- 
resolution ultraviolet, extreme ultraviolet and X-ray tele- 
scopes onboard the SOHO and TRACE satellites have 
revealed the structure of the solar corona as highly fil- 
amentary and inhomogeneous. In particular, TRACE al- 
lowed observations of coronal structures be made down 
to spatial scales smaller than 1000 km and provided ev- 
idence for a fine structuring o f coronal loops consist- 
ing of many individual threa ds ( Aschwanden et al. 2000t 
McEwan fc de Moortell l2006l ). These loop filaments are 
very thin with a thickness of the order of the spatial 
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resolution element of TRACE. Moreover, radio propaga- 
tion measurements have shown that the outer corona is 
also highly inhomogeneous in the direction perpendicular 
to the magnetic field. This nonuniforniity appears in the 
form of filamentary ray-like structures, extending radially 
from the coronal base into the corona. The perpendicu- 
lar length scale of the associated density filaments can 
be as small as 1 k m at the coronal base, and is about 
10 km at 2-5 Rg (|Wool [l99 6. 2006), which is by about 
two orders of magnitude smaller than the observational 
limit of TRACE. Strong inhomogeneity generally prevails 
in the lower corona, and precisely at the boundaries be- 
tween dilute open funnels and dense closed loops, which 
can be associated with sharp gradients in the background 
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plasma quantities. Such gradients play a crucial role in 
the theory of waves in nonuniform plasma. The interac- 
tion of waves with plasma inhomogeneities brings impor- 
tant new physical effects, such as dispersion, phase mix- 
ing (iHevvaerts fc Priesti Il983h and resonant absorption 
(jlonsonlllQsl K 



Additionally, if the plasma has gradients of the density, 
temperature, pressure or magnetic field a plasma current 
may exist and thus provide free energy for driving micro- 
instabilities. They can be very important in the coronal 
context, since they may constitute an affluent source of 
the high-frequency ion-cyclotron waves that have been in- 
voked to play a pro minent role in coronal heating through - 
kinetic dissipation ( Kohl et al.|[l997 : Cranmer et al. I I1999I: 



Marsch fc Tul 120011 ). but for which the c o ronal origin 



still remains unclear. lAxford fc McKenzid (|l992l . ^mB) 



suggested that ion-cyclotron waves could originate in 
the lower corona from small-scale reconnection events 
in the magnetic network. In the extended corona ion- 
cyclotron waves may be locally generated through a tur- 
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or via mode conversion driven by reflection or r efrac- 
tion of low- frequency MHD waves (jMatthaeus et al ][l999 



Cranmer fc van Ballegooijenl 20051 ) . Another possible see 
nario of local generation of ion-cyclotron waves is by 
plasma micro-instabilities that are driven by an inter- 
mittent electron heat flux and plasma outflows accomp a- 
nying micro-flares events jMarkovskii fc Hollwed[2004bH : 
Voitenko fc Goossend 2002 ). by cross-fi eld current fluctu- 
ations of low-frequency MHD modes ( Markovskiil I2OOI 



Markovskn fc Hollwed 1200 2^: 'Zhangl l2003l ). or bv coronal 
ion beams (jMecheri fc Ma rsch 2003). 



The aim of this paper is to study the possible oc- 
currence of micro-instabilities in the ion-cyclotron fre- 
quency range, which are induced by cross-field currents 
and occur under typical coronal-hole conditions. In our 
present case, these currents are generated by an ion- 
electron cross-field drift that is supported by a density 
stratification perpendicular to the ambient magnetic field. 
Linear mode analysis is used in the framework of a coUi- 
sionless multi-fluid model. By neglecting the electron iner- 
tia, this model permits the consideration of ion-cyclotron 
wave effects that are absent from the one-fluid MHD 
model. The ion-cyclotron wavelengths are of the order of 
the ion inertial length, Ai = VAi/^i ~ c/ui, i.e., A ~ A^, 
where Vai = Bq/ ^JJiQUiml is the Alfven speed of the ion 
species i with mass rrii and density n^, and Q,i its cy- 
clotron frequency, respectively, uji plasma frequency, and 
where -So(= \/ -\- B^^) denotes the background mag- 
netic field magnitude. In the weakly coUisional corona, 
the length A^ is much smaller than the electron-ion coUi- 
sional mean free path Ag^, i.e., A^ <C Ag^, and consequently 
A <C Aei, which justifies the coUisionless limit adopted for 
this study. 



Fig. 1. Funnel magnetic field geometry as obtained from 
the 2- D potential field model derived bv lHackenberg et al 



( 2000f ). The field lines emerge from the boundary between 
two adjacent supergranules {x = 0) and expand rapidly 
to fill the corona. The photospheric level is at z = 0. 
Coordinate axes, wave vector and drift velocity are shown 
in red. 

Realistic models of the density and temperature, as 
well as a 2-D funnel model describing an open mag- 
netic field, are used to define the background plasma. 
Considering the WKB approximation (i.e., the wave- 
lengths of interest are smaller than the non-uniformity 
length scale), we first solve locally the dispersion relation 
and then perform a non-local wave analysis using the ray- 
tracing theory. This theory allows us to compute the ray 
path of the unstable waves in the funnel as well as the 
spatial variation of their growth rate. 

This paper is structured as follows. In SectO we 
present the 2-D analytical funnel model used in this study 
to represent an open-field region in a coronal hole. Then 
in Sect [31 we describe how the local and non-local (ray- 
tracing) linear perturbation analysis is carried out, using 
the multi-fluid model. The results are presented and dis- 
cussed in Sectm and finally we give our conclusions in 
SectEl 

2. Background plasma configuration 

For the backgrou nd plasma dens i ty and temperature we 
use the model of iFontenla et aL (Il993 ) for the chromo- 
sphere and the model of Gabriell ~{ 1976 ) for the lower 
corona. The 2-D potential-field and curr ent-free funnel 
model derived bv lHackenberg et"al] ( 2000t ) is used to de- 
fine the background magnetic field (Fig[T]). Analytically, 
this model is given by: 
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arctan- 
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with the relevant parameters given as follow: D=30 Mm, 
d=OM Mm, Soo=11.8 G, Bmax = l-5 kG. 

To model nonuniformity, we consider the length scales 
of the background plasma density and pressure in the 
z-direction, H, to be much bigger than the one in the 
x-direction, L, i.e., H — {{l/nQj)dnoj/dz)~^ ^ L — 
{{l/nQj)dnQj /dx)^^ . Thus only a density gradient in the 
x-direction will be considered, for which use the following 
analytical density model: 



noj{x) = noj{l 



L' 



(3) 



whereby it is assumed that the density increases linearly in 
the x-direction, but only locally within the studied region 
where it is getting bigger by a factor 10. This assumption 
avoids a large growth of the density while the x-coordinate 
increases. 

3. Linear perturbation analysis 

The fluid equations associated with the particle species j 
are given by: 

dn 



dt 



+ V • (rijVj ) = 0, 



(4) 



mi 



9vj 



V, • Vvj) + Vp, - (7,n,(E + V, x B) = 0, (5) 



dp 
'di 



dt 



(6) 



where m-,-, n^, v^, pj and 7j are respectively the mass, 
density, velocity, pressure (which for simplicity is here as- 
sumed to be isotropic) and the polytropic index. The sub- 
script i stands for electron e, proton p or alpha particle a 
(He^"'"). The electric field E and the magnetic field B are 
linked through Faraday's law: 



9B 



3.1. Linearization procedure 



(7) 



We express all the quantities in the above equations as 
a sum of an unperturbed stationary part (with subscript 
0) and a perturbed part (with subscript 1) much smaller 
than the stationary part as follows: 

rij = nQj{x) + nij, pj = poj{x) + , 

= + , B = Bo(a:;, z) + Bi, E = Eq -f Ei, 

nij < noj,pij <_poi, |vij| < Wd]\, |Bi| < |Bo| , (8) 

where we have considered charge neutrality in the unper- 
turbed stationary plasma, i.e., J^j 9j"-0j = 0. Due to den- 
sity stratification in the x-direction, a current is carried by 



drifting electrons and ions in opposite directions parallel 
to the y-axis, with a drift velocity given by: 



1 c% 



L jj^j Bqz' 



(9) 



where Csj = ^J-fjkBTj/rrij 



and — 



q-jBo/rUj 



are, re- 



spectively, the acoustic speed and the cyclotron frequency 
of the species j, and ks is the Boltzman constant. 

The zero-order terms cancel out when Eq.® is in- 
serted into the multi-fluid Eqs.Q-®. Neglecting the 
products of first-order nonlinear terms, we get a system of 
linear equations: 



i(a; — k 
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(12) 



where all the perturbed quantities have been expressed in 
terms of the amplitudes of plane waves. Fourier analysis 
has turned all derivatives into factors: d/dt — > —iuj and 
V — !■ ik, where uj is the wave frequency and k the wave 
vector. Note that in Eq. (11) the term Vij ■ Vv^j has 
been neglected in comparison with the term (k • VDj)'Vij 
since it is smaller by the ratio 2ttLbo,/X where Lbq^ — 
(|Vi3oz|/-Soz)~^ is the length scale of the z-component of 
the background magnetic field and A the wave length of 
interest, i.e.: 

vij • Vvz5j/(k • VDj)vij « L-^lJ {2Tr/\) 

« A/(6Lb„J « 10-^ « 1. (13) 

In this context we want to mention, and remind the reader 
of, the classical guiding-center drift caused by the gravity 
force (g is the gravitational acceleration) in a magnetic 
field, i.e.: 



'93 



1 
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In the corona this drift would lead to a current perpen- 
dicular to the magnetic field, and could thus also provide 
free energy which, as was shown in considerable detail by 
Brinca et al. ( 20021 2003 ). could lead to various kinds of 
plasma instabilities and wave excitation. In comparison 
with Eq. Q , we obtain as an order-of-magnitude estimate 
for the ratio of these two drift speeds the result: 



VDj 



g QjL 



L 



(15) 



where we used again the gravitational scale height H as 
defined in the previous section for the background density 
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being thermally stratified in the z-direction. So the small 
horizontal density variation considered here yields a larger 
drift than the one induced by vertical gravity. However, 
transverse currents driven by gravitational stratification 
may play a role in generating plasma instability at the 
gyrokinetic scale, but the growth rate is expected to be 
smaller by the ratio H/L. 

3.2. Dispersion relation 

The above linearized equations are combined in order to 
obtain a linear relation between the current density Ji 
and the electric field Ei: 



x=9.4 Mm, z=2.5 Mm, 99 85.3°, 9 = 85° 



VDe=fDp=0 



■ Slow mode 
— Alfven mode 
Fas! mode 



Jl — (T • El, 



(16) 



where <r is the conductivity tensor which is related to the 
dielectric tensor e through the following relation: 




kVAp/n,, 



(17) 



The dispersi on relation is obtained using electrodynamics 
theory (e.g.. IStixjlToQ^ : 



D{uj,'k,v) = Dot 



— k X (k X E) + e(cj, k, r) • E 



where c is the speed of light in vacuum and r is the large- 
scale position vector. 

3.3. Ray-tracing equations 



Fig. 2. Two-fluid dispersion curves at the funnel location 
(x=9.4 Mm, z=2.5 Mm) with a Bo-inclination angle ip = 
85.3°. The angle of propagation is 6* = 85°. Top: The case 
of a uniform plasma density, i.e. v^j — 0, with a plasma 
beta f3e = Pp ^ 0.0097. Bottom: In the presence of a drift 
= 0, (18) due to a density gradient with a scale length of L=l km. 

Here w and k are normalized, respectively, to the proton 
cyclotron frequency, fip, and the proton inertial length, 
^p/Vap, where Vap = Bq / ^ ^onopTrip is the proton Alfven 
speed (/io is the magnetic permeability in vacuum). Here 

Toe = Top. 



Considering the WKB approximation, the ray-tracing 
problem consists in solving a system of ordina ry diffcr- 
ential equ ations of the Hamiltonian form ([Weinb erg 1962; 
Lighthil]|ll978h . These equations, that solve for the wave 
frequency u, the wave vect or k and the space coordinat e 
r, have been formulated bv lBernstein fc Friedland ( 1984), 
and in the simple case of a Herniitian dielectric tensor are 
given by: 



dw _ dDiuj,k,r)/dt 
"¥ ' ^ dD{uj,k,r)/dLU 

dk _ dD{uj,-k,r)/dr 
dt ^ dD{uj,k, r)/duj' 

dr _ 

dt ~ ~ 



0, 



dD{uj,k, r)/dk 
' dD{uj,k,r)/duj 



(19) 



(20) 



(21) 



A generalization of these equations to the case of an 
anti-Hermitian dielect r ic ten sor was also proposed by 
Bernstein fc Friedlandl ( 1984 ). In this case, additionally 
to the ray path, the growth rate of the instability can also 
be computed. Note that Eq. (fT9|) is set to zero since the 
dispersion relation does not depend explicitly on the time 
t (the background plasma is stationary) . The above differ- 
ential equations represent a set of initial values problem 
which can be solved using initial conditions obtained from 
the local solutions of the dispersion relation (EqfT8| . 



4. Numerical results 

4.1. Local stability analysis 

The dispersion relation Eq. (|18p is solved numerically for 
the funnel location (x=9.4 Mm, z=2.5 Mm) characterized 
by a magnetic field inclination angle of (p = 85.3° with 
respect to the normal on the solar surface. The value of the 
density length scale in the x-direction, L, is chosen to be 
1 km according to results o btained from radio propagation 
measurements (IWoolll99fil I2OO&) . We consider the cases 
of a two-fluid (e-p) and a three-fluid (e-p-He^"*") model, 
where we consider the effect a second ion population of 
alpha particles (He^+) with typical coronal abundance. 
First, the local solutions of the dispersion relation Eq. ([T5|) 
are presented, and then the results obtained from the non- 
local ray-tracing equations. 



4.1.1. Two-fluid (e-p) drift-plasma 

The two-fluid drift-plasma configuration consists of a 
plasma made of electrons with density riQe and protons 
{riQp) in relative drift to each other in opposite directions 
along the y-axis and perpendicular to the background 
magnetic field with velocities given by Eq. ([§]). The drift 
velocity is sustained by the presence of a local density gra- 
dient in the x-direction, as described by Eq.([3]). For the 
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purpose of comparison, the dispersion curves in the case 
of a vanishing density gradient, L~^=0, and thus no rela- 
tive drift between electrons and ions, i.e., vjji,=vjjp=0, are 
given in the top panel of Figl^l corresponding to an angle 
of propagation 9 = 85°. Consequently, there is no free en- 
ergy in the plasma, and all the modes remain stable. In 
this case the dispersion relation Eq.(fT8| is quadratic and 
three stable modes are present. Each one is represented 
by an oppositely propagating (w > and uj < 0) pair of 
waves. These modes represent the extensions of the usual 
Slow (dotted line), Alfven (dashed line), and Fast (solid 
line) MHD modes into the high-frequency domain around 
a; = f2p {=eBo/mp), where these waves become dispersive. 

However, in the case L=l km, the relative drift be- 
tween electrons and protons is not zero, and the disper- 
sion curves are strongly modified with a breaking of the 
symmetry between forward {uj > 0) and backward {lu < 0) 
propagating waves, as it is shown in the bottom panel of 
FiglH This behavior affects particularly the slow mode. 
Indeed, the free energy provided by the relative drift be- 
tween electrons and protons, which is sustained by the 
density gradient, leads to the appearance of regions of in- 
stability resulting from the coupling between forward and 
backward propagating slow modes. These two single ini- 
tially stable modes couple into one single unstable mode, 
for which the normalized growth rate j/^p is given by the 
red curve, which does not exceed 0.007. 

As it is shown in the top panel of FiglH for small angle 
of propagation, this instability is very weak and concen- 
trated at small wave number kVAp/^^p «C 1, and extends 
to larger k with sensitively bigger growth rate 7 at large 
angle of propagation 9. This instability is also shown in 
the bottom panel of Figl3]as a function of the propagation 
angle 9 and the length scale of the density gradient L in 
the x-direction at the location (x=9.4 Mm, z—2.5 Mm). 
It is clearly seen, that with increasing L, the instability 
growth rate decreases and covers gradually a smaller range 
in the angle of propagation 9. The maximum growth rate is 
7/Op « 0.008 for L=l km and for 9 between 75° and 85°. 
For larger L, the growth rate decreases and the instabil- 
ity appears only at large angles of propagation {9 > 85°). 
This result is expected since the larger the length scale of 
density L is the more decreases the particles drift speed, 
by virtue of (Eql9|) , providing thus less free energy to the 
plasma for driving the instability. 

4.1.2. Three-fluid (e-p-He2+) drift-plasma 

In the three-fluid drift-plasma configuration, in addition 
to electrons (noe) and protons (nop), we consider the pres- 
ence of a second population of ions, namely alpha particles 
He^"'' (indicated by a) with a typical coronal abundance of 
noa — O.lnop. The electrons and ions are in opposite rela- 
tive drift according to Eq.® and caused by the presence 
of a density gradient in the x-direction as in Eq.Q. 

For the purpose of comparison, the dispersion curves 
in the case of L^^~Q (zero gradient of density) and con- 



x=9.4 Mm, z=2.5 Mm, ip w 85.3°, 9 = 85' 



VDe=VDp=VDa 



SJuH'mode 1 




kVAp/n„ 

Fig. 4. Three-fluid dispersion curves at the location 
(x=9.4 Mm, z=2.5 Mm) with a Bo-inclination angle (p = 
85.3°. The angle of propagation is 9 = 85°. Top: The case 
of a uniform plasma, i.e., voj = 0, with the plasma-beta 
f3e = 1.2Pp = 4,13a ~ 0.0097. Bottom: In the presence 
of a drift due to a density gradient with a length scale 



L=l km. Here Tqq — 3T( 



Op 



3Tr 



Oe- 



sequently vue—VDp—VDa—^, are again given in the top 
panel of FigUl corresponding to the location (x=9.4 Mm, 
z=2.5 Mm) with p = 85.3° and to an angle of propaga- 
tion 9 = 85°. In this case the dispersion relation Eg. p^ 
is quadratic and the dispersion diagrams show the pres- 
ence of five stable modes. Each one of them is represented 
by an oppositely propagating (w > and w < 0) pair of 
waves. We note the presence of: slow modes 1 and 2 (re- 
spectively dotted and dashed-dot-dot lines), intermediate 
modes 1 and 2 (dashed and dashed-dot lines), and one 
fast mode which is not shown here because it presents a 
cut-off frequency around ujco ~ O.Gfip, which is out of the 
fre quency range presented in this study (for further details 



Mecheri fc Marschl2006() . Consequently, the considera- 
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tion of a second population of ions leads to the appearance 
of an additional slow and ion-cyclotron mode. 

In the case a non-zero density gradient in the x- 
direction characterized by a length scale L=l km, the dis- 
persion relation is not quadratic anymore, and the sym- 
metry between forward and backward propagating modes 
is broken (bottom panel of FigH]). Indeed, the free energy 
provided by the relative drift between electrons and ions 
sustained by the density gradient leads to the appearance 
of regions of instability. The first instability is similar to 
the one found in the two-fluid model It results from the 
coupling between forward and backward propagating slow 
modes 1. The second instability results from the coupling 
between forward and backward propagating slow modes 2 
with a sensitively smaller growth rate, 7 < 0.001517^ as 
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x=9.4 Mm, z=2.5 Mm, if « 85.3° 





Fig. 3. Growth rate of the two-fluid instability resulting form the coupling between forward and backward propagating 
slow modes at the funnel location (x=9.4 Mm, z=2.5 Mm) with a Bo-inclination angle oiip — 85.3°. Top: As a function 
of the angle of propagation 9 and the normalized wave number, kVAp/^p, for a density length-scale in the x-direction 
of L=\ km {vop = VDe ~ 0.0057Vap). Bottom: As a function of 9 and L for kVAp/^p — 1. 

x=9.4 Mm, z=2.5 Mm, if « 85.3° 




90° { 







2 4 6 


8 




kVAp/Ilp 






1 2 3 4 5 

L(km) 

Fig. 5. Growth rate of the three- fluid slow mode 2 instability which is resulting from the coupling between forward 
and backward propagating slow modes and is absent in the two-fluid case. The results are given for the funnel 
location (x=9.4 Mm, z=2.5 Mm) with a Bo-inclination angle oi tp — 85.3°. Top: As a function of the angle of 
propagation 9 and the normalized wave number kVAp/^p for a density- length scale in the x-direction of L=l km 
(w_De — i-'^VDp — 4ui3q, « 0.0057yAp). Bottom: As a function of 9 and L for kVAp/^p = 1. 
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compared to the first one. This instability appears in gen- 
eral at very large angle of propagation, i.e. 9 > 84° (top 
panel of FigO- As 9 increases, the growth rate appears 
in two different ranges of the wavenuniber k, the first at 
small k, i.e. kVAp/^p < 3, and the second at larger k, i.e. 
kVAp/^p ^ 5. The instability growth rate in general in- 
creases sensitively between kVAp/^p ~ 4 and 8, and has 
its maximum for kVAp/^p ~ 5. However, it is approxi- 
mately 6 times smaller than the slow-mode- 1 maximum 
instability growth rate. When plotted as a function of the 
propagation angle 9 and the length scale of the density 
gradient L (bottom panel of FigE]), this instability mainly 
appears at large angle of propagation, i.e. 9 > 87°, and 
the corresponding growth rate decreases with increasing 
L. The maximum growth rate occurs for L—l km and for 
a propagation angle around 88°. 

4.2. Non-local stability analysis 

In this section we perform a non-local wave study using 
the ray- tracing equations presented in Sect. (|3.3| ). These 
equations are computed using initial conditions obtained 
from the local solutions of the dispersion relation Eq. 
at the starting location (xo=9.4 Mm, yo=0 Mm, zo=2.5 
Mm) . We consider both the two- fluid (e-p) and the three- 
fluid (e-p-He^+) case with a constant concentration of the 
alpha particles (He^"*"), i.e. = O.lrip along the ray path. 
The density length scale in the x-direction L is chosen to 
be 1 km. 

Two-fluid (e-p) drift-plasma 

As presented in the previous section, in the two-fluid case 
the local solution of the dispersion relation showed the 
presence of an instability resulting from the coupling be- 
tween forward and backward propagating slow modes. The 
ray path of this unstable wave, as well as the spatial vari- 
ation of its growth rate 7, when the wave is launched at 
the initial location (xo=9.4 Mm, yo=0, zq~2.5 Mm), are 
presented on Fig[51 The results are shown for different ini- 
tial angle of propagation 9q , with which a different initial 
normalized wave number kg is associated that is chosen 
to always correspond to the maximum growth rate jmax- 
The wave starts propagating from its initial location and 
is well guided along the field lines up to approximately 
the location xwll Mm and z«2.6 Mm, where its motion 
starts to become unguided and irregular (see Fig[6^, [6)d 
and[6t). This unguided motion continues until the wave 
reaches a maximum height where it is reflected and starts 
to propagate downward almost in the (y,z) plane to reach 
again its initial height, i.e., z—2.5 Mm. This maximum 
height reached by the unstable wave is larger for smaller 
initial angle of propagation ^o- Additionally, the smaller 
^0 is the more deeply can the wave propagate in the y- 
direction. The corresponding instability growth rate 7 is 
very irregular, and increases and decreases intermittently 
along the ray path (Fig[6li and[6]3). The growth rate 7 



shows the appearance of peaks of maximum instability, 
which are more important for large initial propagation 
angle 6*0 but do not exceed 7 = 0.02557^. However, the 
instability growth ceases at approximately the maximum 
height reached by the wave, from where on it starts prop- 
agating downward. 

Three-fluid (e-p-He^+) drift-plasma configuration 

In this three-fluid model case, and at the location 
(x=9.4 Mm, z==2.5 Mm), the solution of the dispersion 
relation shows the presence of two instabilities both re- 
sulting from the coupling between forward and backward 
propagating slow modes. For the first instability, which is 
similar to the one obtained in the case of the two-fluid 
model, the ray path as well as the spatial variation of its 
growth rate 7 show a similar behavior as in the two-fluid 
case. The ray path as well as the growth rate variation of 
the second instability, when the wave is launched at the 
initial location (xo=9.4 Mm, yo=0, zo=2.5 Mm), are pre- 
sented in FiglT] The results show that the unstable wave 
starts to propagate as a guided mode for a very small 
distance but then is rapidly reflected to propagate down- 
ward (FigllK, lib andjTh). Additionally, the waves propa- 
gates only to a very short distance in the y-direction, i.e. 
y < 0.004 Mm (FigEt). 

5. Conclusion 

We have studied gradient-drift electromagnetic instabili- 
ties in a coronal funnel, using a coUisionless two- fluid (e-p) 
and three- fluid (e-p-He^+) plasma model with finite pres- 
sure. While neglecting electron inertia, this model allows 
us to take ion-cyclotron wave effects into account. We have 
considered a density stratification transverse to the am- 
bient magnetic field with the typical small length scales 
suggested by the observations. 

First, a local perturbation analysis has been per- 
formed, whereby the dispersion relation has been solved 
for a local region in a coronal funnel. By comparison with 
the results obtained for a uniform-density plasma, it was 
found that the dispersion curves are strongly modified. 
For certain ranges of the wave number, two initially stable 
modes merge into a single unstable mode. Indeed, the free 
energy provided by the relative drift between the different 
plasma species due to the density gradient leads to the ap- 
pearance of regions of instability. In the two-fluid model, 
this instability results from the coupling between forward 
and backward propagating slow modes. For large angles of 
propagation this instability extends over a wide range in 
wave number and is characterized by a small phase speed 
and an increasing growth rate as the angle of propagation 
(with respect to the normal on the solar surface) increases. 

In addition, the consideration of a second ion popula- 
tion of alpha particles (He^+) with typical coronal abun- 
dance led to the appearance of a new instability which also 
results from the coupling of oppositely propagating slow 
mode waves. However, this second instability in general 
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Fig. 6. Ray trajectory and growth rate of the two-fluid slow modes instability. The waves are launched at the initial 
location (xo=9.4 Mm, yo=0, zo~2.5 Mm) in the coronal funnel (with a Bo-inclination angle of ipo « 85.3°) and for 
different initial angle of propagation and normalized wave number ko, chosen for each 9o at the maximum growth 
rate. The dashed lines represent the funnel magnetic field lines, (a) 3-D ray path; (b) Projection of the ray path on 
the x-z plane; (c) Ray path projected on the y-z plane; (d) Growth rate as a function of the z-coordinate; (e) Growth 
rate as a function of the y-coordinate. 



covers only a small range in the wavcnumber domain and 
exists mainly for very large angles of propagation > 82°, 
and it has a substantially smaller growth rate by approx- 
imately a factor 6 as compared to the first instability. 

The non-local analysis, which has been performed 
using ray-tracing equations, revealed that the unstable 
waves, from their launching point, start to propagate up- 
ward for a very small distance but then are refiected and 
propagate downward. During this propagation, the corre- 
sponding instability growth rate has a variable character 
and is increasing, decreasing and vanishing intermittently. 

Consequently, drift currents caused by the fine struc- 
turing of the density in the magnetically open funnels of a 
coronal hole can provide enough energy for driving plasma 
micro-instabilities. They may in turn constitute a prolific 
source of the high-frequency ion-cyclotron waves that have 
been invoked to play, through kinetic wave dissipation, a 
prominent role in heating the open corona. 
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